Dynamics of an inhomogeneous quantum phase transition 



Jacek Dziarmaga^ and Marek M. Rams^ 

^ Institute of Physics and Centre for Complex Systems Research, 
Jagiellonian University, Reymonta 4, 30-059 Krakow, Poland 

We argue that in a second order quantum phase transition driven by an inhomogeneous quench 
density of quasiparticle excitations is suppressed when velocity at which a critical point propagates 
across a system falls below a threshold velocity equal to the Kibble-Zurek correlation length times 
the energy gap at freeze-out divided by h. This general prediction is supported by an analytic 
solution in the quantum Ising chain. Our results suggest, in particular, that adiabatic quantum 
computers can be made more adiabatic when operated in an "inhomogeneous" way. 
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I. INTRODUCTION 

A quantum phase transition is a qualitative change in 
the ground state of a quantum system when one of the 
parameters in its Hamiltonian passes through a critical 
point. In a second order transition a continuous change is 
accompanied by a diverging correlation length and van- 
ishing energy gap. The vanishing gap implies that no 
matter how slowly a system is driven through the transi- 
tion its evolution cannot remain adiabatic near the crit- 
ical point. As a result, after the transition the system 
is excited to a state with a finite correlation length ^ 
whose size shrinks with increasing rate of the transition. 
This scenario, known as Kibble-Zurek (KZ) mechanism 
(KZM), was first described in the context of finite tem- 
perature transitions [H 01 • Although originally motivated 
by cosmology KZM at finite temperature was con- 
firmed by numerical simulations of the time-dependent 
Ginzburg-Landau model Q and successfully tested by 
experiments in liquid crystals [3] , superfluid helium 3 , 
both high- Tc 16,] and low-Tc tZi] superconductors, and con- 
vection cells 0. Most recently, spontaneous appearance 
of vorticity during Bose-Einstein condensation driven by 
evaporative cooling was observed in Ref. However, 
the quantum zero temperature limit, which is in many 
respects qualitatively different, remained unexplored un- 
til recently, see e.g. Refs. |10l-l20l|. The recent interest is 
motivated in part by adiabatic quantum computation or 
adiabatic quantum state preparation, where one would 
like to cross a quantum critical point as adiabatically 
as possible, and in part by condensed matter physics of 
ultracold atoms, where it is easy to manipulate param- 
eters of a Hamiltonian in time and which, unlike their 
solid state physics counterparts, are fairly well isolated 
from their environment. In fact, an instantaneous quench 
to the ferromagnetic phase in a spinor BEG resulted in 
finite-size ferromagnetic domains whose origin was at- 
tributed to KZM [2l|. However, since the transition rate 
was effectively infinite in that experiment, the KZ scal- 
ing relation between the average domain size ^ and the 
quench rate has not been verified. 

The KZM argument is briefly as follows When a 

transition is driven by varing a parameter g in the Hamil- 



tonian across an isolated critical point at gc, then we can 
define a dimensionless distance from the critical point as 
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(1) 



When e — > the correlation length ^ in the ground state 
diverges as ^ ^ kl "i ^ind the energy gap A between 
the ground state and the first excited state vanishes as 
A ~ Setting h= 1 from now on, a diverging A~^ ~ 

|e|~^'' is the shortest time scale on which the ground state 
can adjust adiabatically to varying e(i). A generic e(t) 
can be linearized near the critical point e = as 
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where the coefficient tq is called a quench time. Assum- 
ing that the system was initially prepared in its ground 
state, its adiabatic evolution fails at a e when the time t 
left to crossing the critical point equals the shortest time 
scale A~^ on which the ground state can adjust. Solving 
this equality, we obtain 



'Q 
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(3) 

(4) 



From e the evolution becomes impulse, i.e. the state does 
not evolve but remains frozen in the ground state at e, 
until — e when the evolution becomes adiabatic again. In 
this way, the ground state at e with a KZ correlation 
length 

e - ^ (5) 

becomes the initial excited state for the adiabatic evolu- 
tion after — e. In particular, determines density of 
quasiparticles excited during the phase transition 

d ^ (6) 

in D dimensions. Note that when tq is large, then e is 
small and the linearization in Eq. ([2]) is self-consistent 
because the KZM physics happens very close to the crit- 
ical point between — e and -l-e. 
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II. 



INHOMOGENEOUS TRANSITION 



As pointed out in the finite temperature context |22| . 
in a realistic experiment it is difficult to make e exactly 
homogeneous throughout a system. For instance, in the 
superfluid ^He experiments [5| the transition was caused 
by neutron irradiation of helium 3. Heat released in 
each fusion event, n + ^He — )■ *He, created a bub- 
ble of normal fluid above the superfluid critical temper- 
ature Tc- As a result of quasiparticle diffusion, the bub- 
ble was expanding and cooling with a local temperature 
T{t,r) = exp(-rV2-Dt)/(27rDt)3/2, where r is a distance 
from the center of the bubble and I? is a diffusion coef- 
ficient. Since this T{t,r) is hottest in the center, the 
transition back to the superfiuid phase, driven by an in- 
homogeneous parameter 



e{t,r) 



(7) 



proceeded from the outer to the central part of the bubble 
with a critical front rdt), where e = 0, shrinking with a 
finite velocity v = drddt < 0. 

A similar scenario is likely in the ultracold atom gases 
in magnetic/optical traps. The trapping potential results 
in an inhomogenous density of atoms p{r) and, in general, 
a critical point gc depends on atomic density p. Thus 
even a transition driven by a perfectly uniform g{t) will 
be effectively inhomogeneous. 



eit,r) = 



9c[pir)] 



(8) 



with the surface of critical front, where e = 0, moving 
with a finite velocity. 

According to KZM, in a homogeneous symmetry 
breaking transition, a state after the transition is a mo- 
saic of finite ordered domains of average size ^. Within 
each finite domain orientation of the order parameter is 
constant but uncorrelated to orientations in other do- 
mains. In contrast, in an inhomogeneous symmetry 
breaking transition |22| . the parts of the system that 
cross the critical point earlier may be able to commu- 
nicate their choice of orientation of the order parameter 
to the parts that cross the transition later and bias them 
to make the same choice. Consequently, the final state 
may be correlated at a range longer than ^, or even end 
up being a ground state. In other words, the final den- 
sity of excited quasiparticles may be lower than the KZ 
estimate in Eq. (jS]) or even zero. 

From the point of view of testing KZM, this inhomoge- 
neous scenario, when relevant, may sound like a negative 
result: an imperfect inhomogeneous transition supresses 
KZM. However, from the point of view of adiabatic com- 
putation or adiabatic state preparation it is the KZM 
itself that is a negative result: no matter how slow the 
homogeneous transition is there is a finite density of exci- 
tations © which decays only as a small power of transi- 
tion time tq. From this perspective, the inhomogeneous 



transition may be a way to suppress KZ excitations and 
prepare the desired final ground state adiabatically. 

To estimate when the inhomogenuity may actually be 
relevant, in a similar way as in Eq. ([2]), we linearize the 
parameter e(i, n) in both n and t near the critical front 
where e{t, n) = 0: 



e{t, n) f=a a {n — vt) 



(9) 



Here n is position in space, e.g. lattice site number, a is a 
slope of the quench and v is velocity of the critical front. 
When observed locally at a fixed n, the inhomogeneous 
quench in Eq. Q looks like the homogeneous quench in 
Eq. ^ with 



1 
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(10) 



The part of the system where n < vt, or equivalently 
e(t, n) < 0, is already in the broken symmetry phase. 
The orientation of the order parameter chosen in this 
part can be communicated across the critical point not 
faster than a threshold velocity 



(11) 



When V ^ V the communication is too slow for the inho- 
mogenuity to be relevant, but when v <^ v we can expect 
the final state to be less excited than predicted by KZM. 

Given the relation pO)) . the condition (ITT|) can be 
solved either as 



a 1+" 



(12) 
(13) 



or as a relation between the threshold transition time and 
the slope, 



TQ 



a 1+" 



(14) 



This relation means that, for a given inhomogenuity a, 
the transition is effectively homogeneous when tq ^ tq , 
but the inhomogenuity becomes relevant when the tran- 
sition is slow enough, tq 3> tq. In the homogeneous 
limit a — >■ 0, the threshold transition time fg — ^ oo. 

The threshold velocity in Eq. (|lll) appeared for the 
first time in the context of finite temperature classical 
phase transitions where it looks formally the same, 
but the underlying physics is qualitatively different: the 
scales ^ and t are determined not by the gap of a quantum 
Hamiltonian, but the relaxation time of an open classical 
system. Nevertheless, the key mechanism that non-zero 
order parameter penetrates from the symmetry broken 
phase into the symmetric phase ahead of the critical front 
seems to be the same. 

In the next Section we rederive results (|12ll3ll4p from 
a different perspective. 
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III. 



KZM IN SPACE 



References |23|, l2J] considered a "phase transition 
is space" where e(n) is inhomogeneous but time- 
independent. In the same way as in Eq. this pa- 
rameter can be hnearized in n — Ur, 



e{n) w a (n — 12^) 



(15) 



near the static critical front at rt = where e = 0. The 
system is in the broken symmetry phase where n < Uc 
and in the symmetric phase where n > Uc- In the first 
"local approximation" , we would expect that the order 
parameter behaves as if the system were locally uniform: 
it is nonzero for n < ric only, and tends to zero as (uc — 
n)^ when n — > n~ with the critical exponent /3. However, 
this first approximation is in contradiction with the basic 
fact that the correlation (or healing) length ^ diverges as 
^ \e\~'^ near the critical point and the diverging ^ is the 
shortest length scale on which the order paramater can 
adjust to (or heal with) the changing e(n). Consequently, 
when approaching n~ the local approximation (ric — n)^ 
must break down when a local correlation length ^ ~ 
[a(nc — n)]~'^ equals the distance remaining to the critical 
point {ric — n). Solving this equality with respect to ^, 
we obtain 



a 



(16) 



From n — ric ~ — ^ the evolution of the order parameter in 
n becomes "impulse" , i.e, the order parameter does not 
change until n — ric — +S, in the symmetric phase where it 
begins to follow the local e(n) again and quickly decays to 
zero on the same length scale of 1^. This "KZM in space" 
predicts that a nonzero order parameter penetrates into 
the symmetric phase to a depth 



5n 



a 1+" 



(17) 



The critical point is effectively "rounded off" on the 
length scale of ^. As a consequence, we expect a non- 
zero gap scaling as 



(18) 



as opposed to the local approximation, where we would 
expect gapless quasiparticles near the critical point. 

We expect the finite gap in Eq. to prevent exci- 
tation of the system even when the critical point ric in 
Eq. (|24p moves with a finite velocity, ricit) = vt, up to a 
threshold velocity 



A- 



a 1+" 



(19) 



which is identical with the v in Eq. ([T3|. 

In the following Sections we test these predictions in 
the quantum Ising chain. 
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FIG. 1: The critical front in Eqs. (I31I46|) . 



IV. QUANTUM ISING CHAIN 

The model is 



N 



N-l 



(20) 



For — > cxD, a uniform system with g„ = g has two 
critical points at g = ±1 separating a ferromagnetic 
phase, when \g\ < 1, from two paramagnetic phases, 
when \g\ > 1. We focus on the critical point at g = 1 
when e = g — 1. Given z — I and ly = 1, we expect 



1 



(21) 



independent of either tq or a. The quench is exactly 
solvable for homogeneous g [H, [lH| , but even in an in- 
homogenous case some useful analytic insights can be 
obtained as follows. 

After Jordan- Wigner transformation to spinless 
fermionic operators c„, = 1 — 2cl^Cn and — 
- (cn + 4) n™<„(l " 2ct„c„,), Eq. (HOI) becomes 

N N-1 

H = 2 ^ g„ clcn - ^ + c„+ic„ + h.c.)(22) 

n— 1 n— 1 

This quadratic H is diagonalized to H — J2m ^mlmlm 
by a Bogoliubov transformation c„ — X]m=o('"n"i7m + 
'^nmlm) witli TO numerating N eigcnmodcs of stationary 
Bogoliubov-de Gennes equations 



(23) 



with LOm > 0. Here u^„ 



V. ISING CHAIN: KZM IN SPACE 

To begin with, we consider the ground state of the 
quantum Ising chain in a static inhomogeneous tranverse 
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field gn which can be linearized near the critical point 
g = 1 as 



e(n) 



9r, 



1 



a {n- Uc) , 



(24) 



compare with Eq. ([T5|) . The chain is in the (broken 
symmetry) ferromagnetic phase where n < ric and in 
the (symmetric) paramagnetic phase where n > ric- We 
want to know if the nonzero ferromagnetic magnetization 
Zn = {<^n) the ferromagnetic phase penetrates across 
the critical point into the paramagnetic phase and what 
is the depth Sn of this penetration. 

Since in a homogeneous system quasiparticle spectrum 
is gapless at the critical point only, we expect low energy 
quasiparticle modes ^ to be localized near the critical 
point at He where we can use the linearization in Eq. 
(l24t . We also expect that these low energy modes are 
smooth enough to treat n as continuous and make a long 
wavelength approximation 



d 



n=pl,m 



/"T zr 7/ 



T 

n.m 



(25) 



in Eq. ([25)1 . Under these assumptions, we obtain a long- 
wavelength equation 

WmUm = 2Q!(n - nc)u^ ± 29„u^ . (26) 

After some algebra, its eigenmodes can be found as 



Ura{n) OC tp. 
Vmin) CX Ip. 



l{x) + 1pm{x) , 

i(a;) - -ipmix) , 



where 



— y/a{n — ric) 



(27) 



(28) 



is a rescaled position, 'ipm>aix) are eigenmodes of a har- 
monic oscillator satisfying 



1 



(29) 



and ip-ilx) — 0. As expected, the modes in Eq. (P7)) are 
localized near n = ric where x = 0. A typical width of 
the lowest energy eigenmodes is <5a; ~ 1, or equivalently 



-1/2 



(30) 



When a <C 1 then Sn 1 and the long wavelength 
approximation in Eqs. (|25l26p is self-consistent. Thus 
6n in Eq. (pO)) is the relevant scale of length near Uc and 
we expect that this 5n determines the penetration depth 
of the spontaneous ferromagnetic magnetization into the 
paramagnetic phase. 

We test this prediction by a numerical solution for an 
inhomogeneous transverse magnetic field 



1 



tanh [a{n — ric)] 



(31) 
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FIG. 2: In A and B, exact numerical spontaneous magnetiza- 
tion as a function of n — ric and x = y^(n — Uc) respectively. 
In C and D, ferromagnetic correlation between the site 100 in 
the ferromagnetic phase and a site n when ric = vt — 150 as 
a function of n — vt and in Eq. H36p respectively. Results 
in C and D were obtained with the Vidal algorithm [2^ for 
a = and N = 256. When v ^ 2 there is no ferromagnetic 
correlation across the critical point at n — vt = 0, see C, and 
when i; ^ 2 the correlation penetrates into the paramagnetic 
phase to a depth of Jx„ ~ 1, see D. 



which is shown in Fig. ([T]) with a variable slope a. This 
field can be self-consistenly linearized near n = ric as in 
Eq. (|24p because, when the slope a ^ 1, the predicted 
6n ~ much shorter than the width a ^ of the 

tanh. 

FigureslHA. and B show how the spontaneous ferromag- 
netic magnetization Z„ = (cr^) from the ferromagnetic 
phase, where n < ric, penetrates into the paramagnetic 
phase, where n > ric- In particular, the collapse of the 
rescaled plots in Fig. demonstrates that the penetra- 



-1/2 



as pre- 



tion depth is (5a; ~ 1 equivalent to (5n ~ a 
dieted in Eqs. p7l30p and Ref. [l^. Paramagnetic spins 
near the critical point are biased towards the direction of 
spontaneous magnetization chosen in the ferromagnetic 
phase. 

Moreover, the analytic solution ((?7|) implies a finite 
(relevant) gap 



A 



Wo 



8a 



(32) 



in accordance with the scaling ^ a^/^ predicted by the 



general Eq. p8)) and the numerics in Ref. [23[. This 
gap is the energy of the lowest relevant (even parity) 
excitation of two quasiparticles. 
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VI. ISING CHAIN: INHOMOGENEOUS 
TRANSITION 

Let the critical front in Eq. (|24l) and Fig. [l]niove with 
a velocity v > 0: 

nc{t) = vt . (33) 
A t-dependent version of the long- wavelength Eq. (l26l) , 



idt 



[2a{n - vt)(T^ + 2iayd„ 



,(34) 



can be solved exactly for both v < 2 and v > 2 with 
qualitatively different solutions in the two regimes. Not 
incidentally, v = 2 is the maximal velocity of quasiparti- 
cles at the critical point whose dispersion is uj — 2\k\ for 
small \k\ ^ TT. 



Case of i; < 2 



When V < 2 equation p4l) has solutions 



Vmit,n) cx e-'"'"* 



2 \ 3/4 



a m , 



(35) 



where m — 0,1, 2, the phase ip = arcsin(u/2)/2, and 

„2\ -1/4 



1 



V 

T 



/a (n — vt) 



(36) 



is a rescaled position. When v we recover the static 
solutions (P7)) . In the reference frame of Xv, which is 
co-moving with the critical point, the solutions pSI) are 
stationary modes with > so there are no quasipar- 
ticles in the system, 



d{v <2) = , 



(37) 



and, in particular, no kinks where (;„ = 0. 

As shown in Figs. [2]3 and D, ferromagnetic correla- 
tions penetrate across the critical point into the param- 
agnetic phase to a depth Sxy ~ 1 equivalent to 



2\ 1/4 



^1/2 



(38) 



The penetration depth driy shrinks to when w — >■ 2~ 
suggesting communication problems across the critical 
point when v > 2. 

The same Suy is a typical width of the lowest eigen- 
modes in the spectrum (|35l). As it shrinks to when 
V ^ 2~ , the eigenmodes become inconsistent with the 
long- wavelength approximation in Eq. (1341) . 



B. 



Case of V > 2 



When V > 2 then equation (IMl) can be mapped to a 
homogeneous transition. Indeed, we replace 



1 - 



t - 



n , 



(39) 



introducing local time t measured from the moment the 
critical point passes through n, and simultaneously make 
a transformation 




(40) 



bringing Eq. p4p to a new form 



idf 



-2 — CT^ 



2icj"dn + -9s 

IV 



(41) 

Here — \J\ — + ^ cr^ . Up to an unimportant 
rotation of a Pauli matrix — >■ cr" and the momentum- 
dependent energy shift -^da, the new Eq. (1411) is a ho- 



mogeneous version of the old Eq. (|34|) . but with a longer 

effective quench time tq — tq (l— :^) "^^ > tq. 

Consequently, a quasimomentum representation 
{u^,u~) = {ak,bk) exp(«A:ri — 4ifci/ti)/-\/27r brings the 
homogeneous Eq. (HIT) to the Landau-Zener form: 



ds 



ak 
bk 



-[ ~dk s a"" + a'' 



ak 
bk 



(42) 



where s = fc i is a new time variable and dk = 1/4/c^-rQ 
is a new transition rate. The Landau-Zener formula 
Pk = exp(— 7r/2(5fe) gives excitation probability for a 
quasiparticle k and density of excited quasiparticles is 



d{v > 2) 



^ dk 



(1 



4 
T7I 



,3/4 



2tiJ2¥q 



(43) 



where A ~ 1 is an ultraviolet cut-off. The integral is 
accurate for tq » 1. When v^2 the density 



d{v > 2) 



271^27^ 



^KZM 



(44) 



is the same as the density after a homogeneous quench 
with the same tq, see Ref. [l^, but when u — )■ 2+ then 
d is suppressed below the "homogeneous" density dKZM 
by the factor (1 - 4,/v'^f/^. 



C. Numerical results 

Since the long-wavelength Eq. (|25p does not give 
self-consistent long-wavelength solutions when w — > 2, 
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FIG. 3: Numerical simulations of A'^ = 400 spins. In A 
final density of kinks ^(tq) for different slopes a, and in B 
a rescaled final density a~^^'^ d{v). The solid lines are guide 
to the eye. The fit in panel A shows that when v ^ 2 then 



d 



^1/2 



'Q 



like in the homogeneous KZM, and when v <^ 2 



then d is suppressed below the homogeneous KZM density. 



velocity of quasiparticles at the critical point. Excitation 
of kinks is dramatically suppressed when a critical front 
propagates slower than v and the ferromagnetic phase is 
able to communicate its choice of ferromagnetic polar- 
ization to the paramagnetic phase ahead of the front. In 
contrast, when the front is much faster than v the com- 
munication across the front is not efficient enough and 
kinks are excited as in the homogeneous KZM. However, 
even above v density of excited kinks is suppressed below 



0.06 



0.05 



0.04 



-0 0.03 



0.02 



0.01 



a = 2 



+ numerics 
i 

KZM 



-d = d_.(l-4/v ) 



2,3/4 



KZM 



OW-14 



10 



we simulated the exact time-dependent version of the 
Bogoliubov-de Gennes equations 



dui 



dt 



2"^Ti.™ ' (45) 



on a finite lattice of N sites for a time-dependent trans- 
verse magnetic field 



9n{t) 



1 



tanh [a{n — vt)] 



(46) 



with a moving critical point at ric = ''^t, compare Eq. (j3ip 
and Fig. [TJ Results are shown in Figs. [3] and HI Figure S] 
demonstrates good quantitative agreement between Eq. 

and numerical results, despite the breakdown of the 
long- wavelength approximation near v = 2. 



FIG. 4: Comparison between Eq. (|43|) (solid blue), the ho- 
mogeneous KZM (dotted red), and numerical simulations on 
a lattice oi N = 1000 spins (crosses) at a fixed slope a = 2~®. 



/ .2\3/4 

the "homogeneous" KZM density by a factor (1 — ^1 

which is significantly less than 1 when v is close to . 

Thus the general estimates (jllll2ll3ll4l) are confirmed 
by our solution of the quantum Ising chain, but we leave 
their interesting implications when z 7^ 1 or in more than 
one dimension for future exploration. 

The estimates and the solution suggest that "inhomo- 
geneous" adiabatic quantum computers can be more adi- 
abatic than their "homogeneous" counterparts. 



VII. 



CONCLUSION 
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